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Simulation of relativistic shocks and associated radiation 
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Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we 
investigated long-term particle acceleration associated with a relativistic electron-positron jet prop- 
agating in an unmagnetized ambient electron-positron plasma. We have also performed simulations 
with electron-ion jets. The simulations were performed using a much longer simulation system than 
our previous simulations in order to investigate the full nonlinear stage of the Weibel instability for 
electron-positron jets and its particle acceleration mechanism. Cold jet electrons are thermalized 
and ambient electrons are accelerated in the resulting shocks for both cases. Acceleration of ambient 
electrons leads to a maximum ambient electron density three times larger than the original value 
for pair plasmas. Behind the bow shock in the jet shock strong electromagnetic fields are generated. 
These fields may lead to time dependent afterglow emission. We calculated radiation from electrons 
propagating in a uniform parallel magnetic field to verify the technique. We also used the new 
technique to calculate emission from electrons based on simulations with a small system with two 
different cases for Lorentz factors (15 and 100). We obtained spectra which are consistent with those 
generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent 
magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel 
instability. 



I. RPIC SIMULATIONS 



Relativistic Jets Injected into Unmagnetized 
Plasmas using a Large System 
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Particle-in-cell (PIC) simulations can shed light on 
the physical mechanism of particle acceleration that 
occurs in the complicated dynamics within relativis- 
tic shocks. Recent PIC simulations of relativistic 
electron-ion and electron-positron jets injected into an 
ambient plasma show that acceleration occurs within 
the downstream jet |H1H 0, H, S M 
ITT1 . [l2l [l3l ]. In general, these simulations have con- 
firmed that relativistic jets excite the Weibel insta- 
bility, which generates current filaments and associ- 
ated magnetic fields [H, [l5[ , and accelerates electrons 



Therefore, the investigation of radiation result- 
ing from accelerated particles (mainly electrons and 
positrons) in turbulent magnetic fields is essential for 
understanding radiation mechanisms and their ob- 
servable spectral properties. In this report we present 
a new numerical method to obtain spectra from par- 
ticles self-consistently traced in our PIC simulations. 



We have performed simulations using a system with 
{L x ,L y ,L z ) = (4005A, 131A, 131A) (A = 1: grid 
size) and a total of ~ 1 billion particles (12 particles 
/cell/species for the ambient plasma) in the active grid 
zones [l6j |. We have performed two kind of simula- 
tions: an electron-positron jet is injected into an ambi- 
ent pari plasma and an electron-ion jet into electron- 
ion ambient plasma (mi/m e = 20). In the simula- 
tions the electron skin depth, A co = c/u pc — 10. OA, 
where ui pc = (47re 2 n e /m e ) 1/ ' 2 is the electron plasma 
frequency and the electron Debye length A c is half of 
the grid size. Here the computational domain is six 
times longer than in our previous simulations 0, H[ • 
The electron number density of the jet is 0.676n c , 
where n e is the ambient electron density and 7 = 15 
for both cases. The electron/positron thermal veloc- 
ity of the jet is v? th = 0.014c, the ion thermal velocity 
where c = 1 is the speed of ligh t. The ion thermal 
velocity of jet is v- th = «? th * yjm~~jm/ l = 0.00313c. 
The electron/positron thermal velocity in the ambient 
plasma is w° th = 0.05c {v\ th = 0.0112c). 

Figure 1 shows the averaged (in the y—z plane) elec- 
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FIG. 1: The averaged values of electron density (a, c) and 
field energy (b, d) along the x at t = 3750WP,, 1 . The left 
panel (Figs, la and lb) is for the case of pair plasmas and 
the right panel (Figs, lc and Id) is for the electron-ion case 
(mi/m c = 20). Figs, la and lc show jet electrons (red), 
ambient electrons (blue), and the total electron density 
(black). Figs, lb and Id show electric field energy (red) 
and magnetic field energy (blue) divided by the total ki- 
netic energy. 



tron density and electromagnetic field energy along 
the jet at 3750WP,, 1 for electron-positron (Figs, la and 
lb) and electron-ion (Figs, lc and Id) jets. The re- 
sulting profiles of jet (red), ambient (blue), and total 
(black) electron density are shown in Figs, la and lc. 
The ambient electrons are accelerated by the jet elec- 
trons and pile up towards the front part of jet. For 



electron-positron jet, at the earlier time the ambient 
plasma density increases linearly behind the jet front. 
At the later time the ambient plasma shows a rapid 
increase to a plateau behind the jet front, with ad- 
ditional increase to a higher plateau farther behind 
the jet front. The jet density remains approximately 
constant except near the jet front. 

The acceleration of ambient electrons becomes vis- 
ible when jet electrons pass about x/A = 500. The 
maximum density of accelerated ambient electrons is 
attained at t — 1750u;^" e . The maximum density 
gradually reaches a plateau as seen in Fig. la. The 
maximum electromagnetic field energy is located at 
x/A = 1, 700 as shown in Fig. lb. 

The Weibel instability remains excited by contin- 
uously injected jet particles and the electromagnetic 
fields are maintained at a high level, about four times 
that seen in a previous, much shorter grid simulation 
system (L x = 640 A). At the earlier simulation time 
a large electromagnetic structure is generated and ac- 
celerates the ambient plasma. As shown in Fig. lb, 
at the later simulation time the strong magnetic field 
extends up to x/A = 2,000. These strong fields be- 
come very small beyond x/A = 2000 in the shocked 
ambient region 

In the case with electron-ion jet, due to the heavier 
ions at 3750cj po 1 the density is piled up only slightly 
(Fig. lc). The generated electromagnetic fields are 
smaller than those for the electron-positron case. Fur- 
thermore, in the front of electron-ion jet the electro- 
magnetic fields disappear as shown in Fig. Id. In 
order to generate a shock it will require at least sev- 
eral time longer simulations. 



II. THE STANDARD SYNCHROTRON 
RADIATION MODEL 

A synchrotron shock model is widely adopted to de- 
scribe the radiation mechanism in the external shock 
thought to be res pon sible for observed broad-band 
GRB afterglows [H M, 0J, HI [S|. Associated 
with this model are three major assumptions that are 
adopted in almost all current GRB afterglow models. 
Firstly, electrons are assumed to be "Fermi" acceler- 
ated at the relativistic shocks and to have a power-law 
distribution with a power-law index p upon accelera- 
tion, i.e. N(E c )dE c oc E~ p dE c . This is consistent 
with recent PIC simulations of the shock formation 
and particle acceleration [l2| and also some Monte 
Carlo models [H [H, S3, but see H Hi]. Sec- 
ondly, a fraction e G (generally taken to be < 1) of 
the total electrons associated with ISM baryons are 
accelerated, and the total electron energy is a frac- 
tion e e of the total internal energy in the shocked re- 
gion. Thirdly, the strength of the magnetic fields in 
the shocked region is unknown, but its energy density 
(B 2 /8tt) is assumed to be a fraction es of the internal 
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energy. These assumed "micro-physics" parameters, 
p, e c and eg, whose values are obtained from spectral 
fits [13, Ha| reflect a lack of knowledge of the under- 
lying microphysics [29]. 

The typical observed emission frequency from an 
electron with (comoving) energy 7 e m e c 2 in a frame 
with a bulk Lorentz factor r is v = Tj 2 (eB /27rm e c). 
Three critical frequencies are defined by three charac- 
teristic electron energies. These are v m (the injection 
frequency), v c (the cooling frequency), and vyi (the 
maximum synchrotron frequency) . In our simulations 
of GRB afterglows, there is one additional relevant 
frequency, z/ a , due to synchrotron self-absorption at 
lower frequencies [2(J Hi], [3(| HH . 

The general agreement between the blast wave dy- 
namics and the direct measurements of the fireball 
size arg ue for the validity of this model's dynamics 
[20l |2l| . The shock is most likely collisionless, i.e. 
mediated by plasma instabilities (29[. The electro- 
magnetic instabilities mediating the afterglow shock 
are expected to generate magnetic fields. Afterglow 
radiation was therefore predicted to result from syn- 
chrotron emission of shock accelerated electrons [32|] • 
The observed spectrum of afterglow radiation is in- 
deed remarkably consistent with synchrotron emission 
of electrons accelerated to a power-law distribution, 
providing support for the standard afterglow model 
based on synchrotron emission of shock accelerated 
electrons [l3 M M M, M M, M, M ■ 

In order to determine the luminosity and spectrum 
of synchrotron radiation, the strength of the magnetic 
field (ee) and the energy distribution of the electrons 
(p) must be determined. Due to the lack of a first 
principles theory of collisionless shocks, a purely phe- 
nomenological approach to the model of afterglow ra- 
diation was ascribed without investigating in detail 
the processes responsible for particle acceleration and 
magnetic field generation [29j. Rather, one simply 
assumes that a fraction £b of the post-shock ther- 
mal energy density is carried by the magnetic field, 
that a fraction e is carried by electrons, and that 
the energy distribution of the electrons is a power- 
law, dlogn e /(iloge = p (above some minimum en- 
ergy £0 which is determined by e e and p), 6b, e and 
p are treated as free parameters, determined by ob- 
servations. It is important to clarify here that the 
constraints implied on these parameters by the obser- 
vations are independent of any assumptions regarding 
the nature of the afterglow shock and the processes 
responsible for particle acceleration or magnetic field 
generation. Any model should satisfy these observa- 
tional constraints. 

The properties of synchrotron (or "jitter") emis- 
sion from relativistic shocks will be determined by the 
magnetic field strength and structure and the electron 
energy distribution behind the shock. The character- 
istics of jitter radiation may be important to under- 
standing the complex time evolution and/or spectral 



structure in gamma-ray bursts [37| . For example, 
jitter radiation has been proposed as a means to ex- 
plain GRB spectra below the peak frequency that are 
harder than the "line of death" spectral index asso- 
ciated with synchrotron emission [38l . [39| . i.e., the 
observed spectral power scales as F v cx ^ 2 / 3 , whereas 
synchrotron spectra are F u tx v 1 ^ or softer [39[. 
Thus, it is essential to calculate radiation production 
by tracing electrons (positrons) in self-consistently 
treated small-scale electromagnetic fields. 



III. NEW NUMERICAL METHOD FOR 
CALCULATING SYNCHROTRON EMISSION 

Let a particle be at position ro(i) at time t pol. l4ll 
l42j . At the same time, we observe the electric field 
from the particle from position r. However, because 
of the finite velocity of light, we observe the particle at 
an earlier position ro (t ) where it was at the retarded 
timei' =t-St =<-R(t')/c HereR(t') = |r-r (t')| 
is the distance from the charge (at the retarded time 
t ) to the observer. 

After some calculation and simplifying assumptions 
the total energy W radiated per unit solid angle per 
unit frequency from a charged particle moving with 
instantaneous velocity (3 under acceleration can be 
expressed as [H, Ell 



d-*W 
dilduj 



16tt 3 



(1) 



n x [(n - /?) x 0] iu (t'- a . ro ( t 'vc),,> 
(1-0 -n)* 



Here, n = R(t )/|R(t )| is a unit vector that points 
from the particle's retarded position towards the ob- 
server. 

The observer's viewing angle is set by the choice of 



n (nl+n; 



1). The choice of unit vector n along 



the direction of propagation of the jet (hereafter taken 
to be the x-axis) corresponds to head-on emission. For 
any other choice of n (e.g., 6 1 = 1/7), off-axis emission 
is seen by the observer. 



A. Synchrotron radiation from two electrons 
propagating in parallel uniform magnetic field 

In order to calculate radiation from relativistic jets 
propagating along the x direction [40j we consider 
a test case which includes a parallel magnetic field 
{B x ), and jet velocity of Vji.2 = 0.99c. Two elec- 
trons are injected with different perpendicular veloc- 
ities (fj_i = 0.1c, Uj_2 = 0.12c). A maximum Lorenz 
factor of 7max = {(1 - (vf 2 + vj_ 2 ) / 'c 2 }- 1 ' 2 = 13.48 is 
calculated with the larger perpendicular velocity. 
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FIG. 2: The case with a strong magnetic field (B x = 3.7) 
and larger perpendicular velocity = 0.1c, v±2 = 

0.12c). The paths of two electrons moving helically along 
the x— direction in a homogenous magnetic field shown 
in the x — y-plane (a). The two electrons radiate a time 
dependent electric field. An observer situated at great 
distance along the n-vector sees the retarded electric field 
from the moving electrons at the rest frame (b). The ob- 
served power spectrum at different viewing angles from 
the two electrons (c). Frequency is in units of uip e . 



Figure 2 shows electron trajectories in the x — y 
plane (red: v±2 = 0.12c, blue: v±% = 0.1c) (a: 
left panel), the radiation (retarded) electric field (b: 
middle panel), and spectra (right panel) for the case 
B x = 3.70. The two electrons are propagating left to 
right with gyration in the y—z plane (not shown). The 



gyroradius is about 0.44A for the electron with the 
larger perpendicular velocity. The seven curves show 
the power spectrum at viewing angles of 0° (red), 10° 
(orange), 20° (yellow), 30° (moss green), 45° (green), 
70° (light blue), and 90° (blue). The higher frequen- 
cies become stronger at the 10° viewing angle. The 
critical angle for off-axis radiation 7 = 180°/(7r7 ma x) 
for this case is 4.25°. As shown in this panel, the spec- 
trum at a larger viewing angle (> 20°) has smaller 
amplitude. 

Since the jet plasma has a large velocity x- 
component in the simulation frame, the radiation 
from the particles (electrons and positrons) is strongly 
beamed along the x-axis (jitter radiation) (3R l39j. 

Equations 6.30a and 6.30b show that the radiation 
with the viewing angle a = disappears (see Fig. 
6.5 in the textbook of Rybicki and Lightman [43[). 
However, based on other textbooks, radiation at the 
viewing angle 0° should not vanish [H, H(| . This 
aspect is shown in Fig. 2c, and at the higher frequency 
the amplitude at the viewing angle 10° is stronger 
than that with viewing angle 0°. 



B. Calculating Synchrotron and Jitter Emission 
from Electron Trajectories in Self-consistently 
Generated Magnetic Field 

In order to validate our numerical method we 
performed simulations using a small system with 
(L x ,L y ,L z ) = (645A, 131A, 131A) (A = 1: grid 
size) and a total of ~ 0.5 billion particles (12 
particles/cell/species for the ambient plasma) in the 
active grid zones Q- First we performed simulations 
without calculating radiation up to t = 450w^ e 1 . The 
jet front is located around about x/A = 480. We se- 
lected 12,150 electrons for each jet and ambient elec- 
trons randomly. Recently, a similar calculation has 
been carried out for the radiation from accelerated 
electrons in laser-wakefield acceleration [47| and in 
shocks [Hj]. 

Figure 3 shows (a) the current filaments generated 
by the Weibel instability and (b) the phase space of 
x/ A—'yVx for jet electrons (red) and ambient electrons 



(blue) at t = 450W" 1 for 7 = 15. 

Figure 4 shows (a) the x-component of current den- 
sity generated by the Weibel instability and (b) the 
phase space of jet electrons and ambient electrons af- 
ter t s = 75u)~* (at t = 525<j~£). 

We calculated the emission from 12,150 electrons 
during the sampling time t s = t% — t\ — 75w~ e 
with Nyquist frequency o->n = 1/2 At = 200cj pe where 
At = 0.0050;"^ is the simulation time step and the 
frequency resolution Auj = l/t s = 0.0133o> pe . 

The spectra shown in Fig. 5 are obtained for emis- 
sion from jet electrons and ambient electrons sepa- 
rately for two case with 7 = 15andl00. In this case 
the spectra are calculated for head-on radiation (0°) 
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FIG. 3: Two-dimensional images in the x — z plane at 
y/A = 65 for t = 450u I 7 e 1 for the case with 7 = 15. 
The colors indicate the x-component of current density 
generated by the Weibel instability, with the x- and z- 
components of magnetic field represented by arrows (a). 
Phase space distributions as a function of x/ A — 7^ x plot- 
ted for the jet (red) and ambient (blue) electrons at the 
same time. 
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FIG. 5: Spectra obtained from jet and ambient electrons 
for the two viewing angles for 7 = 15 (upper) and 7 = 100 
(lower). Spectra with jet electrons are shown in red (0°) 
and orange (5°). Spectra from ambient electrons show the 
lowest levels by blue (0°) and light blue (5°). 
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FIG. 4: Two-dimensional images in the x — z plane at 
y/A = 65 for t = 525<^p e for the case with 7 = 15. 
The colors indicate the x-component of current density 
generated by the Weibel instability, with the x- and z- 
components of magnetic field represented by arrows (a). 
Phase space distributions as a function of x/ A — 7^ x plot- 
ted for the jet (red) and ambient (blue) electrons at the 
same time. 



and 5°. It is noted that In the case with 7 = 100 
spectrum is extended in the higher frequency. How- 
ever, the spectrum with viewing angle 5° is greatly 



decreased in particular in the high frequency due to 
the narrow beaming angle. 

The radiation from jet electrons show 
Brcmsstrahlung-like spectra as a red line (0°) 
and orange line (5°) [41| . The spectra with jet 
electrons are different from the spectra shown in Fig. 
2c. Since the magnetic fields generated by the Weibel 
instability are rather weak and the jet electrons are 
not much accelerated, the trajectories of jet electrons 
are almost straight with only a slight bent. 

We compare these spectra with our known spec- 
tra obtained from two (jet) electrons, the case with a 
parallel magnetic field (B x — 0.37), and jet velocity of 
v ji.2 = 0.99c. Two electrons are injected with different 
perpendicular velocities (v±i — 0.01c, v±2 = 0.012c). 
A maximum Lorenz factor of j max = {(1 — (v? 2 + 
^ 2 )/c 2 } -1 / 2 = 7.114 accompanies the larger perpen- 
dicular velocity. The critical angle for off-axis radia- 
tion Q 1 = 180°/(7T7 max ) for this case is 8.05°. 

Comparing the spectra with Figs. 5 and 6c we 
find similarities. The lower frequencies have flat spec- 
tra and the higher frequencies decrease monotonically. 
The slope in Fig. 5 is less steep than that in Fig. 
6c. This is due to the fact that the spread of Lorenz 
factors of jet electrons is substantial and the average 
Lorenz factor is larger as well. As shown in Fig. 7.16 
in Hededal's Ph. D. thesis [41| . the turbulent mag- 
netic field with the red noise (/x = —3) makes the 
spectrum shifted toward higher frequencies. This ef- 
fect is found in Fig. 5 

We obtained spectra using several different param- 
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FIG. 6: The case with a weak magnetic field (£> x = 0.37) 
and small perpendicular velocity («xi = 0.01c, v±2 = 
0.012c). The paths of two electrons moving helically along 
the x— direction in a homogenous magnetic field shown 
in the x — y-plane (a). The two electrons radiate a time 
dependent electric field. An observer situated at great 
distance along the n-vector sees the retarded electric field 
from the moving electrons (b). The observed power spec- 
trum at different viewing angles from the two electrons 
(c). Frequency is in units of ur- 



eters with jet electrons and ambient magnetic field. 
However, the strength of the magnetic fields gener- 
ated by the Weibel instability is small in the region 
x/A < 500 (cb < 0.07) as shown in Fig. lb, therefore 
the spectra for these cases are very similar to Fig. 5. 



As shown in Fig. 7.12 in Hededal's Ph. D. thesis [411 ] . 
the trajectories of jet electrons have to be chaotic to 
produce a jitter-like spectrum as shown in Fig. 7.22. 

In order to obtain the spectrum of synchrotron (jit- 
ter) emission, we consider an ensemble of electrons 
selected in the region where the Weibel instability has 
fully grown and electrons are accelerated in the gen- 
erated magnetic fields as shown in Fig. 1, which is 
being investigted. 



IV. DISCUSSIONS 

Emission obtained with the method described above 
is obtained self-consistently, and automatically ac- 
counts for magnetic field structures on small scales 
responsible for jitter emission. By performing such 
calculations for simulations with different parameters, 
we can investigate and compare the different regimes 
of jitter- and synchrotron- type emission [381.1391]. The 
feasibility of this approach has already been demon- 
strated [il|, [42j , and its implementation is straight- 
forward. Thus, we should be able to address the low 
frequency GRB spectral index violation of the syn- 
chrotron spectrum line of death 39]. 

Medvedev and Spitkovsky recently showed that 
electrons may cool efficiently at or near the shock 
jump and are capable of emitting a large fraction of 
the shock energy [49|]. Such shocks are well- resolved 
in existing PIC simulations; therefore, the microscopic 
structure can be studied in detail. Since most of the 
emission in such shocks would originate from the vicin- 
ity of the shock, the spectral power of the emitted 
radiation can be directly obtained from finite-length 
simulations and compared with observational data. 

As shown in Fig. 1, behind the trailing shock the 
electrons are accelerated and strong magnetic fields 
are generated. Therefore, this region seems to produce 
the emission that is observed by satellites. We will 
calculate more spectra based on our RPIC simulations 
and compare in detail with Fermi data. 
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